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Abstract: In systems considered for quantum computing, i.e., for control of quantum 
dynamics with the goal of processing information coherently, decoherence and deviation 
from pure quantum states, are the main obstacles to error correction. At low tempera- 
tures, usually assumed in quantum computing designs, some of the accepted approaches 
to evaluation of relaxation mechanisms break down. We develop a new formalism for es- 
timation of decoherence at short times, appropriate for evaluation of quantum computing 
architectures. 
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Consider a microscopic quantum system described by the Hamiltonian H$. This sys- 
tem, S, can be a quantum bit (qubit), or several qubits, controlled externally individually 
and by switching on and off pairwise qubit interactions. This qubit paradigm helps to de- 
fine the questions in describing how S interacts with the surrounding macroscopic world. 
Interactions with the surroundings depend on the setting. For example, in quantum mea- 
surement, the wavefunction of the system is probed, and part of the process involves a 
strong interaction with the measuring device. However, in most applications, the external 
interactions are quite weak. In quantum computing , the aim is to minimize their effect. 

Traditionally, interactions with the surroundings have been modelled by the modes of 
a "bath" B, with the mode Hamiltonians Mjf, 

H b = Y j Mk- (1) 

K 

The interaction of the bath modes with S, will be modelled by 

H i = A s Pb = A s J2 J k, (2) 

K 

where As is a Hermitean operator of S, coupled to the operator Pb of B. For a bosonic 
bath [1,2], which we use as a prototype, 

M K = u) K a} K a K , J K = g* K a K + g K ^ K , (3) 

where the ground state is shifted to zero, and we work in units h = 1. The total Hamilto- 
nian is then H = H$ + H B + Hi. 

Let p(t) represent the reduced density matrix of S, after B has been traced over. For 
large times, the effect of the environment on an "idle" quantum system, i.e., one that is 
not otherwise externally controlled, is expected to be thermalization: this density matrix 
should approach exp (—/3Hs) / Tr s [exp (—/3Hs)], where j3 = 1/kT. For times t > 0, we 
can consider the degree to which the system has departed from coherent pure-quantum- 
state evolution, due to the interactions and entanglement with the bath. The temperature 
and other external parameters needed to characterize p(t), are determined by the properties 
of B, which in turn might interact with "the rest of the universe." 

In the eigenbasis of Hs, we consider the matrix elements of 

H s \n) = E n \n) , p mn (t) = (m\p(t)\n) . (4) 
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For large times, we expect the diagonal p nn to approach values oc e~ l3En , while the off- 
diagonal elements p m ^ n —> 0, corresponding to thermalization and decoherence in the 
energy basis. To establish this, several assumptions are traditionally made [1-5]. At time 
t = 0, it is assumed that the bath is thermalized, i.e., its modes, K, have density matrices 
9k = e~ f3Ml< / Tr k (e~ l3MK ). The overall density matrix R is then unentangled, 

R(t) = P (t)l[e K ] t=o . (5) 



K 



Now, a series of assumptions are made, e.g., the Markovian and secular approxima- 
tions. The Markovian approximation essentially assumes that the density matrices of the 
bath modes are reset to the thermal ones, on time scales shorter than any dynamical times 
of the system interacting with the bath. This amounts to using (5) for times t > and 
is a natural assumption, because each bath mode is coupled only weakly to the system, 
whereas it is "monitored" by the rest of the universe and kept at temperature T. Such 
approaches aim at master equations for p mn at large times, consistent with the Golden 
Rule, thermalization and decoherence. In variants of these formalisms, two external (to 
S) times scales are identified. One is the inverse of the upper cutoff frequency of the bath 
modes, 1/u>d- Another is the "thermal" time K/kT = (3. 

There is evidence [3,5,6] that the Markovian- type approximations are only valid for 
times large than both these time scales. This is hardly a limitation at room temperatures. 
However, for quantum computing, in semiconductor-heterostructure architectures [7-10], 
T ?s 10 — 100 pK. The thermal time then becomes dangerously close to the single-qubit 
control time even for slower qubits, based on nuclear spins. We emphasize that not all the 
approximation schemes have this problem [5] . We also point out that quantum computing 
architectures utilize [7-10] qubits and modes that couple them, that have large spectral 
gaps. It is believed that, especially at low temperatures, spectral gaps slow down relaxation 
processes. Thus, qubit levels are considered in quantum dots, in atoms, in large magnetic 
fields, and coupled by highly nondissipative "quantum" media [8,10], e.g., the quantum- 
Hall two-dimensional electron gas. 

At low temperatures, spectral gaps lead to separation of time scales of the initial deco- 
herence vs. later-stage thermalization and further decoherence driven by energy-conserving 
processes. One can then question whether the energy basis is the best to describe decoher- 
ence for times before the thermalizing processes take over. The issue of the basis has also 
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come up in models of quantum measurement [11,12], where the eigenbasis of the interac- 
tion, As, may be more appropriate. In fact, in quantum computing applications we really 
want to retain a pure quantum state [12,13]. The value of 1 — Tr g [p 2 (t)], of other mea- 
sure, may be more indicative than off-diagonal matrix elements. Therefore, it is desirable 
to have basis-independent expressions for the reduced density operator p(t). 

Recently, several groups have reported [6,10,12,14-16] results for spin decoherence in 
quantum computing systems. Some have not invoked all the traditional approximations, 
Markovian and secular, etc., or have utilized the spectral gap of the bath modes, to achieve 
better reliability of the short-time results. Experimental efforts are picking up momentum, 
with the first limited results available [17] by NMR/ESR techniques. An approach, termed 
adiabatic decoherence, has been developed [12] to avoid the ambiguity of the basis selection 
and achieve exact solvability [6,12,14,15]. The price paid was the assumption that H$ is 
conserved (a version of the quantum nondemolition assumption), i.e., [Hs,H] = [H$, Ag] = 
0. This makes the eigenbasis of H$ and Ag the same, but precludes energy exchange, 
leaving only relaxation pathways that contribute to decoherence. Generally, A5I7) = A 7 |7), 
where the Greek index labels the eigenstates of Ag; the Roman indices will be used for the 
energy, (4), and, capitalized, for the bath modes, (2,3). 

The most widely used approximation has been the second-order perturbative expan- 
sion in the interaction, Hj. In the present work, we report a different scheme, valid for 
short times. It has several advantages, such as becoming exact in the adiabatic case, yield- 
ing several explicit results, and permitting derivation of higher-order approximations. We 
do assume that at t = the system and bath modes are not entangled. Our formulation 
also relies on that the Hamiltonians are all time-independent, and, therefore, applies to 
"idling" (possibly interacting) qubits. It is reasonable to assume that a lower limit on 
decoherence rate can be evaluated in such an idling state. The t = factorization assump- 
tion, shared by all the recent spin-decoherence studies, represents the expectation that 
control by short-duration but large externally applied potentials, will "reset" the qubits, 
disentangling them from the environment to which they are only weakly coupled. Thus, it 
is the qubit system that gets approximately reset and disentangled from the bath at t = 0, 
rather than the bath is thermalized by the "rest of the universe," as assumed in Markovian 
approximation schemes. 

The overall density matrix is R(t) = e - i( - Hs+HB+Hl ^R(0) e i( - Hs+HB+Hl ^. The fol- 
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lowing relation for the exponential factors will be used as our short-time approximation, 

e i(H s +H B +H I )t+0(t 3 ) _ e iH s t/2 e i(H B +Hi)t e iH s t/2 _ /g 



It becomes exact for the adiabatic-decoherence case. Furthermore, e are replaced by 
three consecutive time-evolution-type transformations on R(0). Therefore, the approx- 
imate expression for R(t) will have all the desired positivity properties. Extensions to 
higher-orders in powers of t are possible: see [18] for various expressions valid to 0(t 4 ) and 
0(t 5 ) in our nomenclature. In the resulting approximation to the matrix element, 

p mn (t) = ^BHe-^^e-^^^e-^'/^lOje^^e^^^e^^ln) , (7) 

we insert the decomposition of the unit operator in the S space in terms of the eigenbasis 
of As before the second exponential, and in terms of the eigenbasis of H$ after it, etc., 

Pmn(t)= J2 ltB[e- 4 ^'/ 2 H7><7b>e- 4(Hfl+ ^ ft,) 'e- < ^*/ 2 p M (0) 

7 p q S 



^Y[e K y E ^ 2 e^ HB+x ^\q\5)(8\n)e iE ^ 2 ] . (8) 

K 

We use (1,2) to write 

p mn (t)= Yl {^ Eq+En - Ep - Em)t/ \Mi){i\p)mmppM 



jpq5 



JJ Tr K ~ e -i(M K +\^j K )t Q Ke i(M K +\ s j K )t^ j _ 



X 

K 

This expression actually allows rather straightforward calculations in typical quantum- 
computing applications which involve single or few two-state systems. The overlap brackets 
are between the eigenstates of Hs (labeled by m, n, p and q) and of Ag (labeled by 7 and 
6). The traces are over the modes of B. Since these modes have identical structure, with 
independent coupling constants, the calculation needs only be done once, in the space of 
one mode. 

For the bosonic bath [2], see (3), with the thermal initial 6 K = (l - e~ l3uK ) e "^ a W 5 
the product of traces in (9) is known [6,12,14], 

Pmn(t)= £ {e J ^ + ^-^- E -^ 2 (m|7)(7|p)(?|5)(5|n)p pq (0) 

7 p q S 



|2 



x exp ( — ^ ® K 2 2 (A 7 — Xs) 2 sin 2 coth ^ K + z (A 2 — A 2 ) (sinu^t — cu^t) 
^ ^ u k L 2 2 

(10) 

The last term in the exponent, linear in t, is the "renormalization" of the system energy 
levels due to its interaction with the bath modes. It could be removed by adding the term 
Hr = A| J2k \9k\ 2 /^k to the total Hamiltonian. The non-negative real spectral sums, 
B(t) and C(t), 

B *( t ) = 8 £!^sin^coth^, (11) 



K 



K 



I 12 

C(t) = V (u;** - sinu K t) , (12) 

K K 

when converted to integrals over the bath mode frequencies, with the cutoff at ujd, have 
been discussed extensively in the literature [2,6,14], for several choices of the bath mode 
density of states and coupling strength g as functions of the mode frequency. In summary, 
we get the approximation 



Pmn{t) = £ je^ 
JpqS \ 



-\-E n —E p — E m )t/2 



H7><7b><«IW*|n>ft*(0) 



x exp 



- A B 2 (t) (A 7 -A 5 ) 2 + iC{t) (A 2 -A 2 ) 



(13) 



which is exact for all times in the adiabatic case [6,12,14,15,19], and has the properties of 
a density operator. 

A basis-independent representation for p(t) is obtained by using ^pa exp[— £? 2 (AA) 2 /4] = 
f^dy e~ y2 exp[iy£?(AA)]. Exponential factors in (13) can then be reproduced by applying 
operators on the wavefunctions, and the sums carried out, 

^ p = Jdy e -V 2 e - iH st/2 e i{yB(t)A s +C(t)A 2 s ] e ~iH s t/2 ^ e iH s t/2 e -i[yB(t)A s +C(t)A 2 s ] e iH s t/2 _ 

(14) 

Decoherence is explicit in (14): if p(0) is a projection operator |^o)(^o| (a pure state), 
then p(t > 0) is obviously a mixture (integral over y) of projectors \ip(y, t)){ip(y,t)\, where 

_ e -iH s t/2 e i[yB(t)A s +C(t)A 2 s ] e ~iH s t/2 

As an application, let us consider the case of H$ proportional to the Pauli matrix a z , 
e.g., spin-1/2 in magnetic field, and As = cr x , with the proportionality constant absorbed 
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in the couplings gx in (3). We study the deviation of the state of a spin-1/2 qubit, initially 
in the energy eigenstate | | ) or | |), from pure state, by calculating Tr 5 [p 2 (t)]. For a 
two-by-two density matrix, this trace can vary from 1 for pure quantum states to 1/2 for 
maximally mixed states. With p(0) = |T)(tl or |l)(l|> 

Tr s [p 2 (t)] = l[l + e-^] . (15) 

As the time increases, the function B 2 {t) grows monotonically from zero [2,6,12,14]. For 
Ohmic dissipation, B 2 (t) increases quadratically for short times t < 0{1/ujd), then loga- 
rithmically for 0{l/u D ) < t < 0(h/kT), and linearly for t > 0(h/kT). (For other bath 
models, it need not diverge to infinity at large times.) This calculation illustrates that the 
present approximation can yield reasonable results for short and even intermediate times. 

In summary, we have derived short-time approximations, (13,14), for the density ma- 
trix or its energy-basis matrix elements, for the case of the bosonic heat bath with initially 
fully thermalized modes. Other baths can be studied by using (9). 
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